arXiv:l509.05180v2 [cond-mat.str-el] 17 Nov 2015 


Magnetic impurity in a Weyl semimetal 


Jin-Hua Sun/’^ Dong-Hui Xu,^ Fu-Chun Zhang,and Yi Zhou^’^ 

^Department of Physics, Zhejiang University, Hangzhou 310027, China 
^Collaborative Innovation Center of Advanced Microstructures, Nanjing 210093, China 
^Department of Physics, Hong Kong University of Science and Technology, Clearwater Bay, Kowloon, Hong Kong, China 

(Dated: November 18, 2015) 

We utilize the variational method to study the Kondo screening of a spin-1/2 magnetic impurity 
in a three-dimensional (3D) Weyl semimetal with two Weyl nodes along the fc^-axis. The model 
reduces to a 3D Dirac semimetal when the separation of the two Weyl nodes vanishes. When the 
chemical potential lies at the nodal point, /r = 0, the impurity spin is screened only if the coupling 
between the impurity and the conduction electron exceeds a critical value. For finite but small /r, the 
impurity spin is weakly bound due to the low density of state, which is proportional to , contrary 
to that in a 2D Dirac metal such as graphene and 2D helical metal where the density of states is 
proportional to |/i|. The spin-spin correlation function Juu(r) between the spin n-component of the 
magnetic impurity at the origin and the spin u-component of a conduction electron at spatial point 
r, is found to be strongly anisotropic due to the spin-orbit coupling, and it decays in the power-law. 

The main difference of the Kondo screening in 3D Weyl semimetals and in Dirac semimetals is in 
the spin x- (y-) component of the correlation function in the spatial direction of the z-axis. 

PACS numbers: 75.20.Hr, 03.65.Vf, 71.27.-|-a 


I. INTRODUCTION 


Two-dimensional (2D) Dirac fermions have been pro¬ 
posed and observed in graphene^^^ and in surface 
states of 3D topological insulators (TIs)'^^®. 3D Dirac 
semimetal represents a new state of quantum matter, 
which can host 3D Dirac fermions in the bulk^. The 
stable 3D Dirac semimetals have been realized experi¬ 
mentally in NaaBi compounds® and CdaAsa crystals®’^^, 
where the Dirac points are stabilized by crystalline sym¬ 
metry, and the Dirac nodes are degenerate. If the in¬ 
version (V) or time-reversal (T) symmetry is broken, 
each Dirac node splits into two Weyl nodes resulting 
in Weyl semimetals^^^^®. Weyl semimetals show inter¬ 
esting physics such as Fermi arc surface states and chi¬ 
ral anomaly. Recently, Weyl semimetals have attracted 
much attention because a new TaAs family of Weyl 
semimetals has been predicted in theories^"*’’^® and sub¬ 
sequently observed in experiments^®^^®. 

A single magnetic impurity^^ in a conventional metal 
is well described by the Anderson impurity model. This 
model or the Kondo problem^® has been widely studied 
by using various methods^®^®^. The impurity spin-1/2 is 
fully screened, and the correlation between the impurity 
spin and a conduction electron of distance r is of a power 
law decay l/r‘^ ii r < and l/'r‘^+^ if r > ^k, with 
the Kondo coherence length and d the dimensionality 
of the host metal®®^®^. Graphene has a peculiar electron 
structure, where the density of states (DOS) vanishes at 
a charge neutral point or at half-hlling. The property of 
a magnetic impurity in half-filled graphene falls into the 
category of pseudo-gap Kondo problem®®^''’® which has 
been studied constantly using various methods including 
numerical renormalization group (NRG). In single layer 
graphene, the full screening of a magnetic impurity re¬ 
quires a hnite strength of the hybridization between the 
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FIG. 1. (Color online). Schematics of the dispersion relation 
of (a) Dirac semimetals and (b) Weyl semimetals along the 
fcz-axis. The Dirac cones are located at kz = 0 in the Dirac 
semimetals, and the Weyl nodes are located at kz = ±6/(2A) 
in the Weyl semimetals, p is the chemical potential and eu is 
the energy level of the Anderson impurity. 


impurity and Dirac electrons, and the spin-spin correla¬ 
tion between the impurity and conduction electron de¬ 
cays with 1/r® power law for large Recently the 

Kondo effect in 3D Dirac and Weyl systems was studied 
using the NRG method, and it was found that the mag¬ 
netic impurity shows a diverse range of Kondo physics, 
depending on the DOS of the host system and the sym¬ 
metries broken by perturbations^^. The spin-spin corre¬ 
lation in the spin-orbit coupled systems may be interest¬ 
ing since the coupling between the spin and the momen¬ 
tum results in anisotropy in both the spin and spatial 
spaces. 

The purpose of this paper is to investigate the proper- 
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ties of spin-spin correlation in Dirac and Weyl semimet¬ 
als. We study the binding energy and the various com¬ 
ponents of the spin-spin correlation and we illustrate 
the similarity or differences between the Dirac and Weyl 
semimetals. The variational method we apply has been 
used to study the ground state of the Kondo problem 
in conventional metals^^’^^, antiferromagnets'^^'’ and 2D 
helical metals^^. 

The paper is organized as follows. We present the 
model and dispersion relation in Sec. II. In Sec. Ill, 
we apply the variational method to study the binding 
energy. In Sec. IV, we investigate the spin-spin correla¬ 
tion between the magnetic impurity and the conduction 
electrons in Dirac or Weyl semimetals, and we compare 
the results in these two systems. Finally, the discussions 
and conclusions are given in Sec. V. 

II. ANDERSON MODEL HAMILTONIAN 


two Weyl nodes due to the broken T or parity symmetry. 
We note that the Weyl semimetal studied in the present 
paper corresponds to the T-broken case. For 6 > 0, the 
Weyl nodes are located at the points (0,0, ±6/2A) on the 
/cz-axis as shown in Fig. 1 (b). 

The local magnetic impurity Hamiltonian can be writ¬ 
ten as 

Hd = [sd - + d^di) + Undxndi- (3) 

and are the creation and annihilation operators 
of the spin-up (spin-down) state on the impurity site, 
and Ed and U are the impurity energy level and on-site 
Coulomb repulsion, respectively. 

Finally, the hybridization between the magnetic impu¬ 
rity and the host material is described by 

Hv = Y.div^^. ( 4 ) 

k 


We utilize the Anderson impurity model to study the 
Kondo screening of a spin-1/2 magnetic impurity in a 
3D Dirac or Weyl semimetal. The model Hamiltonian 
contains three parts: the kinetic energy Hq of the Dirac 
or Weyl semimetal, the impurity Hamiltonian Hd, and 
the hybridization between the impurity and the Dirac or 
Weyl semimetal Hy- The Hamiltonian reads 


dk = (dkfjdki)'^, where dks (d^s) is the impurity anni¬ 
hilation (creation) operators in the momentum space. V 
is the hybridization strength and we assume the mag¬ 
netic impurity is equally coupled to the four bands in the 
semimetal. 


{ Vk 0 Vfc 0 \ 
Vo Vk 0 Vk )■ 


( 5 ) 


H = Ho + Hd + Hv. 
Here the kinetic energy part is given by 


( 1 ) The single particle eigen-energy, ej(k), with j = 
1 , 2 ,3,4 can be obtained by diagonalizing the non¬ 
interacting Hamiltonian Hq, and it is given by 


i^o = X] ^k[2Aa^(S X k) • £ + 2Xzaykz + 
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( 2 ) 

where k = {kx,ky,kz) and k± = kx ^ iky. a, S are 
Pauli matrices in orbital and spin spaces, respectively. 
The basis vectors of the bulk states are chosen as = 
(okt, ak 4 ,, ^)kt) ^ki)^j where a and b are orbital indices. 
Mk = M — tk^, with M the Dirac mass. For > 0 
and M < 0, Hq describes a strong topological insulators 
with Z 2 index (1; 000). Here we shall consider the case 
with M = 0 while the bulk gap closes at k = 0 and a 
transition between a topologically nontrivial phase and 
a trivial phase occurs. For simplicity, we shall assume 
that Xz = X eliminates the extra anisotropy, and we set 
Mk = 0 since the nonzero Mk merely modifies the high 
energy dispersion which has a minor influence on our 
study. At 5 = 0, the Hq describes the Dirac semimetal. 
At 6 ^ 0, it describes a Weyl semimetal with broken 
time-reversal (T) symmetry.®’^^’^® In a Dirac semimetal, 
the Dirac points are doubly degenerate [see Fig. 1 (a)], 
while in a Weyl semimetal each Dirac point splits into 


£i(k) = -^4A2(fc2 + fc 2 ) + {b- 2XkzY - At, 


£2(k) ~ ~\J 4A^(fc^ -l- ky) -|- (6 -l- 2Xkz)'^ — fJ., 

! - ( 6 ) 

e3(k) = y' 4A2(fc2 + k^) + {h- 2XkzY - 

e 4 (k) = y^ 4 A 2 (/c 2 - 1 - k^) -I- (6 -I- 2XkzY - 


where 6 = 0 and 6 7 ^ 0 correspond to the Dirac and Weyl 
semimetals, respectively. The corresponding eigenstates 
are given by 


1 


7kl = 


Vcl 

1 


{(/ikakt + 2iAfc_aki + i/’ik 6 kt + 2Afc_6k;}, 


7k2 = ^={-((>2kakt - 2iAfc_ak; -b i(/2k6kt + 2Afc_6k4,}, 
VD 2 

7k3 = ^^{(/>3kakt + 2iAfc_aki + i/'3k6kt + 2Afc_6k;}, 

vD3 

7k4 = ^^{-(()4kakt - 2iAfc_ak; -b i/'4k6kt + 2Afc_6k4,}, 
vD4 

( 7 ) 

where C/ are the normalization factors. The parameters 
^i)/k {I = 1, 2, 3,4) are given by 


(/>ik = (6 - 2Xkz) - TOk, />2k = (6 -b 2Xkz) - Uk, 
</ 3 k = (^ - 2 Afc^) -b TOk, /'4k = (6 -b 2Xk,z) + n-i^, 
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where mk(?^k) = + ky) + (& T ‘2\kz). We can 

rewrite the total Hamiltonian H in the diagonal basis as 

^ - /^)7kj7kj + ^k(7kjt^ki + H.C.) 


ki 


kj 


+ {sd - m) X! ^ 


where the impurity operators dkj are given by 

j (1 + i)(/'ik j , (1 + i)2Afc_ j 
“ki —-7=—Mt 4-7^=- 






4,= 


'Jc^ 


'JC^ 


(1 + i)^()3k j , (1 + i)2Afc_ 

«k3 = -7=-Mt 4-7^=-“ii 




VC~s 


-(l-i)(/'4kj , (l-i)2Afc_ 

“k4 = - 7 = -«t 4-7^=-“i- 






(9) 


( 10 ) 


III. THE BINDING ENERGY 


We start from the simplest limit in which the magnetic 
impurity and the host material are decoupled from each 
other, namely Hy = 0. The ground state wavefunction 
of Hq can be written as 
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FIG. 2. (Color online). Calculated binding energy At of a 
magnetic impurity in Dirac or Weyl semimetal as a function 
of r for three values of chemical potential p. A is the energy 
cut-off, and F = 3V)?/A^ is the effective hybridization. At 
H = 0, there is a threshold F > Fc = |cii|/A = 7.5 x 10“' for 
the bound state. For p 7^ 0, A(, is finite, but too small to be 
seen in the figure for small F. 



The variational principle requires 


i'i^)o= n 7k/io), 

kl;E3(k)<M 


( 11 ) 


where the product runs over all the states inside the 
Fermi surface, li Sd < < Sd + U, the impurity is singly 

occupied with a local moment, and the impurity energy 
is Ed — The total energy of the system is given by 
the sum of the energies of the host material and of the 
magnetic impurity, 

Eo = ed-fJ.+ ^ ej(k)- (-^ 2 ) 

kj;E3(k)</r 


In the above equations, the index j = 1,2, 3,4 denotes 
the four bands in the Dirac or Weyl semimetal. 

Now we study the effect of Hy which describes the hy¬ 
bridization between the magnetic impurity and the host 
Weyl semimetal by using a variational method. We shall 
assume large-17 limit. In this case, we may neglect the 
doubly occupied impurity states in our trial wavefunction 
for the ground state, which is given by 

1 ^) = ico + Yl Ckj4j7k4)|0). (^3) 

kj 

Co and Ckj are variational parameters to be determined 
by optimizing the ground state energy. The energy of the 
Hamiltonian H in the trial state Ifh) is given by 

Ekj [(^0 - £3 (k) + n)clj + 2I4coCkj + (£3 (k) - n)cl] 

c§ + Ek,< 


dE/dcf) = dEjdc]ij = 0, (15) 

from which we obtain two equations below: 

Eco = ^[VkCkj + (ej (k) - /i)co], 

kj (16) 

(^0 “ 1 “ 

which lead to 

(E - Xl(^f (*^) “ ^kCkj) 

kj kj (17) 

[E- Eo + (ej(k) - n)]ck 3 = T4co. 

We define the binding energy as At = Eq — E. If At > 0, 
then the hybridized state has lower energy than the bare 
state, so that the hybridized state is more stable. From 
Eq. 17 we obtain 


_ _ 14 cq _ 

(ej (k) - m) - Af,’ 


and 


{E - Eo-\- {sd - m))co = ^ VkCkj- 


We then obtain the self-consistent equation 
(£d - m) - Ab = ^ ■ 






(£j(k) - n)-At' 


(18) 


(19) 


(14) 


( 20 ) 

















By solving Eq. 20 numerically, we obtain the binding 
energy A^. From Eq. 18 we can calculate Ckj for given 
k and j. In the calculations, we introduce an energy 
cutoff A, and hence the truncation of momentum kc = 
A/(2A). The summation over momentum in Eq. 20 is 
then replaced by an integration ^ I 

One can see that the binding energy is independent of 
the host system to be the Dirac or Weyl semimetal . 
Although the Weyl nodes are located at = ±5/2A in 
the Weyl semimetal, the dispersion relation around each 
Weyl node is exactly the same as those around the Dirac 
cones, thus the summation over k on the right-hand side 
of Eq. 20 shall also remain the same. We define an 
effective hybridization E = . From Eq. 20 we obtain 


(£d ~ k) ~ + r[ 


A A 

— op + k + Ab 

~ A ^ Ab 


( 21 ) 

Then in the limit of small T and TA < p, — + ST^ we 

have 


. . , u — Ed — TA-I-2rn., , , , , 

AfcRiAexp{-2rv7A-(22) 


If /r 0, the hybridization shall always lead to a finite 
binding energy A{, > 0. However, if p, = 0, Eq. 21 gives 
rise to 


or A 

ed + TA = Ab-—Alln—. (23) 

A Afe 

In the limit of small F, we obtain Ab Ri -|- TA. The 
density of states in the Dirac or Weyl semimetal van¬ 
ishes, hence the binding energy A^ has a positive solu¬ 
tion only when the effective hybridization T is above a 
critical value, F > ■^. This is similar to the case of 
graphene^^AS other 2D helical metals^^. In the con¬ 
text of Kondo physics, it is so called pseudo-gap Kondo 
problem. 

The self-consistent solution of the binding energy is 
plotted in Fig. 2 as a function of the effective hybridiza¬ 
tion F. We fix the value of /r — = 7.5 x 10“^A, and the 

chemical potential is chosen as ^ = 0 and k = ±0.1A. 
While ^ = 0 the system is at half-filling, and we can 
see that the binding energy A;, is non-zero only if the 
effective hybridization is greater than the critical value 
be = kdl/A = 7.5 X 10“®. While ^ 0 the binding 

energy always has a positive value, which can also be 
seen from the analytical results shown in Eq. 22. The 
DOS D{E) oc is much smaller than that in graphene 
near the charge neutral point, so that the screening ef¬ 
fect and thus the binding energy are much smaller than 
those in the graphene case. The asymmetry between the 
k = ±0.1A cases is due to the asymmetry of the impurity 
state between being empty and doubly occupied. 




FIG. 3. (Color online). Spin-spin correlation between mag¬ 
netic impurity and conduction electron along the i-axis (panel 
a), and along the y-axis (panel b). The results are the same 
for the Dirac or Weyl semimetal. The inset in panel (a) illus¬ 
trates spin of magnetic impurity at the origin r = 0 and 
conduction electron spin 5“ at a distance r along the x-axis. 
The parameters are k = —O.OIA, Vk = 0.05A and A = 0.05A, 
and the energy cut-off A is large enough that the value of A 
shall not affect the low-energy physics. All the other spin-spin 
correlations not shown here are zero. 


IV. SPIN-SPIN CORRELATION BETWEEN 
MAGNETIC IMPURITY AND CONDUCTION 
ELECTRONS 

Now we study the spin-spin correlation between the 
magnetic impurity Sd = ^d^erd and the conduction elec¬ 
tron spin Sc = |cf(7c in the Dirac or Weyl semimetal. 
The impurity-spin-conduction-electron-spin correlation 
function is evaluated for p, y^ 0 while the binding energy 
Ab always has a positive value, i.e., there exists Kondo 
screening. We assume that the impurity is located at 
the origin r = 0, and we only consider the simplest case 
in which the spin-spin correlation function is evaluated 
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FIG. 4. (Color online). Spin-spin correlation between a mag¬ 
netic impurity at the origin and a condnction electron at a 
distance r along the 2 :-axis for Dirac semimetal {b = 0) and 
Weyl semimetal of three values of b. The two Weyl nodes 
are at b/X along the fc^-axis. Correlation for spins along z, 
Jzz{r), is independent of b and Jzz(_r)\h = Jxx\b^o = Jyy\b=o, 
while Jxxir) and Jyy{r) are fe-dependent. The parameters are 
H = -O.OIA, Vk = 0.05A and A = 0.05A. 


along the three axes: x-, y- and z-axes. The spin-spin 
correlation function between the magnetic impurity and 
the conduction electron is given by 

J„.(r) = (5“(r)5,"(0)), (24) 


ken due to the displacement of the Dirac cones along 
the /cz-axis. Both the systems have rotational sym¬ 
metry along the z-axis, so we have J„„(r) = Ju'v'i^') 
if u' = TZz{(3)u, v' = 'TZz{l3)v, r' = 'R,z{l3)r where 
TZz{(3) is a rotational operator along the z-axis. Actu¬ 
ally one can demonstrate that both the systems are un¬ 
changed under a combined mirror reflection and time- 
reversal transformation. We may denote the mirror re¬ 
flection with respect to the y — z plane as then 

we have TM^''{x,y,z) = {-x,y,z), {kx,ky,kz) = 

ikx,-ky,-kz) eindrMy^iSx,Sy,Sz) = i-Sx,Sy,Sz). 
Thus if r is on the x-, y- or z-axis, we can easily find 
that only Jyz{ry) = -Jzy{ry) and Jxz{rx) = -Jzx{rx) 
are nonzero among the off-diagonal terms. We should 
emphasize that due to the absence of the spin SU(2) sym¬ 
metry, the off-diagonal terms can generally be non-zero 
in the coordinate space, but we may only concentrate on 
the three spatial axes in this present paper. 

We can define functions below and use them to simplify 
the spin-spin correlation function in the coordinate space. 


k 

e“*‘'‘’^2A(-ifc+)cfcj 

k 

r,(r)=^e-*^'-^2A(ife_)c,„ 


(25) 






,,^4X^{kl + kl) 




~ <2^3 ^3- 


where u,v = x,y,z and (...) denotes the ground 
state average. We start with simple symmetry analy¬ 
sis. As we have mentioned, the Dirac semimetals we 
study in this paper preserve both V and T symme¬ 
tries, while in Weyl semimetals the T symmetry is bro- 


Here again j = 1, 2, 3,4 are used to denote the four bands 
in the Dirac or Weyl semimetals. Then the diagonal 
terms and the non-zero off-diagonal terms of the spin- 
spin correlation function along the three axes are given 
by 


■7zz(r) — —{[|di + dap -|- \l2 + / 4 P + \Qi + QaP + 1^2 + >/4p — | Ji -I- J 3 P — | J 2 -I- J 4 P — |Ti -I- Tap — IT 2 -I- T 4 P]}, 
Jxxir) = -{2[(Ji + Ja)(r* + Tl) + (J 2 + Ji)iT* + T*) + {h + h){Q*2 + Ql) + ih + h){Ql + Q^)]}, 

Jyy{v) = -{2[-(Ji + Ja)(T; - Tl) + (J 2 + Ji){n + t;) + (/i + h){Ql + Ql) + {h + h){Q\ + Q^)]}; 

Jxz{'^) = — 2 Re{[/i + + i(/2 + Ii)][Ti +T^ — i{T2 + T^)] — [Ji + J3 — i( J2 + J4:)\{Qi + Qa + i(Q2 + Ql)]}; 

Jy,{v) = -2Im{[/i + /a + i (/2 + /4)][ri* + T; - + Tl)] - [J^ + J3- i(J 2 + Ji)\[Ql +QI+ 1(^2 + Q:)]}- 

(26) 


If we consider the simplest case when y, < 0 and b = 
0 then we have m = 2Xk and = —2A(/c -I- kz) = 
—2Xk{l + cos9), and then the correlation writes 


We take the first term for example 

k 

= E 


ikr^lA: \ '' _ —ikr ^Ik 


-dm(j)ik 


^kl 


-ikr 


Cfei = i^e ‘‘"'■(1-bcos6l)cfei. 


-4mfc 


( 28 ) 


J.. cx I/1P + IJ2P + IQ1P + IQ2P. 


( 27 ) 












After a straightforward integration, we obtain 


1 / ; A I- 1^*1 r-N 

h oc - - e 

r 




+ Af, ;Oi+Ai> /■’■A/2A g-iy 

oT / 


(29) 


2A 


lrAb/2\ V 


Ay 


By assuming fj, + Ai, <C A, we can ignore the second term 
and obtain 

_ 3iT4ao 1 

' 4A (M' ’ 


where oq is the number defined in Eq. 13, and A is the 
energy cut-off. We can easily find that I 2 = Ii, and Qi 
and <52 also decay with l/{kcr)‘^, which indicates that 
the spin-spin correlation Jzz{r) follows a l/{kcr)‘^ decay 
at long distance. Using similar analysis, we can see that 
while 6 = 0 all the other diagonal and off-diagonal com¬ 
ponents of the spin-spin correlation decay with l/(kcr)^ 
in the real space. The results for the Weyl semimetal 
{b ^ 0) are more complicated, which will be discussed 
later. 

In Fig. 3, we show the results of J„„(r) for r along 
the X- and y-axis. The values of Juv{,r) along these two 
axes are independent of b, hence they are the same for 
the Dirac and Weyl semimetals. As shown in Fig. 3, 
the diagonal terms are all antiferromagnetic at short dis¬ 
tance. In Fig. 3 (a) for r along the i-axis, we have 
Jxx{r) = Jyy{r) ^ Jzzir) and only two of the off-diagonal 
terms, Jxzij") = —Jzxir), are non-zero. The inset in Fig. 
3 (a) shows the schematics of the displacement of the 
magnetic impurity and the conduction electrons. We as¬ 
sume the location of the impurity is at r = 0 and the 
conduction electron is on the a;-axis, and the displace¬ 
ment between them is r. and S'" are used to denote 
the spins on the magnetic impurity and the conduction 
electron, respectively. Similarly in Fig. 3 (b) we have 
Jyy — Jzz Jxx and only one off-diagonal term Jyz is 
non-zero. The off-diagonal terms reflect the spin-orbit 
coupling in the Dirac and Weyl semimetals. 

In Fig. 4, we show the spin-spin correlation with r 
along the z-axis, which depends on b, hence is different 
between the Dirac and Weyl semimetals. Firstly, in the 
3D Dirac system, the diagonal terms are all equivalent, 
Jxx = Jyy = Jzz- As the value of b increases, the Dirac 
cones split and Weyl nodes emerge. One can see that 
Jzzir) remains the same as b increases. We always have 
Jzzir)\b=o = Jzz{r)\b>o = Jxx{r)\b=o = Jyy{r)\b=o if r 
is along the z-axis. We still have Jxx = Jyy for 6 > 0, 
and we find that these two terms are modified in the 
coordinate space as b increases. If /i < 0 and r = rz, 
then from Eq. 26 we obtain 


Jxxirz) = -2{JiT* + J 2 T* + hQ* + I 2 QI). (31) 


where the indices 1 and 2 denote the lower bands whose 
Weyl nodes are located at = —6/2A and kz = b/2X, 
respectively. Ji corresponds to the contribution from the 



FIG. 5. (Color online). The product {kcr)^Jxx as a function 
of the dimensionless distance k^r. Here the displacement r 
is along the z-axis. We use p = —O.OIA, I 4 = 0.05A and 
A = 0.05A. The magnetic impurity is coupled to (a) Dirac 
semimetals {b = 0) and (b) Weyl semimetals (6 = 1), respec¬ 
tively. 


lower band whose Weyl node is at kz = —b/2X. If we per¬ 
form a simple translation along the kz axis and substitute 
kz with -I- 6/2A, we obtain Ji = where Jf* 

is the results for Dirac semimetal. Similarly, we can see 
that T* = such that JiT* = jf Tf 

and J 2 T* = Given that the two lower bands 

are degenerate at 6 = 0, we obtain 

Jxxirz) = jJ^(rz)cos(6r/2A), (32) 

where J^^irz) is the x — x spin-spin correlation along 
the z-axis in the Dirac semimetal, and b/2X is half of 
the separation of the two Weyl nodes along the fcz-axis. 
One can easily find that the displacement of the Weyl 
nodes along the fcj-axis will add an extra phase factor 
to each function defined in Eq. 25, and it will further 
induce complexity to the oscillation behavior of the spin- 
spin correlations along the z-axis. One can use a similar 
approach to prove that for the y — y correlation it shall be 
Jyyirz) = Jyyirz) cos(6r/2A). Generally in the short dis¬ 
tance limit, as shown in Fig. 4, the spin-spin correlation 
decays faster as b increases. 

Shown in Fig. 5 is the product ikcr)^Jxxir) as a 
function of the dimensionless distance kcr. We use 
/i = —O.OIA, 14 = 0.05A and A = 0.05A. We can see 
that both the terms for 6 = 0 and b = I.O oscillate in 
coordinate space, and the decay rate is proportional to 
ikcr)~^. This decay rate is consistent with the general 
results in normal metals that the spin-spin correlation de¬ 
cays with id = 3) at far displacement. By carefully 
examining Eq. 26, we find that the effect of b is to add 
an extra phase factor to the spin-spin correlations, i.e., as 
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the phase factor cos(6r/2A) added to the spin-spin corre¬ 
lation J^^{rz) given in Eq. 32. The results for different 
values of b generally have similar contributions, i.e., that 
the oscillation in the coordinate space is modified by b. 

V. DISCUSSION AND CONCLUSION 

In conclusion, we have studied spin-1/2 Anderson im¬ 
purity in a Dirac or Weyl semimetal. We apply the varia¬ 
tional method to study the problem at the large-t/ limit. 
Due to the spin-orbit coupling in these two systems, the 
spatial spin-spin correlations between the magnetic im¬ 
purity and the conduction electron are highly anisotropic. 
The results of the binding energy are quite similar to 
those in graphene and 2D helical metals, and they are 
found to be the same for the Dirac and Weyl semimetal. 
Due to the vanishing DOS at half-filling, there exists 
a critical value of the effective hybridization strength. 
We may obtain a positive binding energy only if the hy¬ 
bridization is above the critical value. While the sys¬ 
tem is particle-hole asymmetric, the DOS at the Fermi 
surface becomes finite, such that the hybridization al¬ 
ways leads to a positive binding energy. However, due 
to the fact that the DOS near the Weyl points is pro¬ 
portional to E^, the screening effect of the conduction 
electrons is much weaker than that in graphene, in which 
DOS{E) oc \E\. Therefore, in the weak hybridization 


limit and at the chemical potential near the Dirac point, 
the impurity binding energy in the 3D Dirac or Weyl 
semimetal is much smaller than that in the correspond¬ 
ing graphene system or 2D helical metal. 

The spin-spin correlations in both the Dirac and Weyl 
semimetal are highly asymmetric in coordinate space due 
to the spin orbit coupling. Although the Kondo temper¬ 
ature of a Dirac or Weyl semimetal is mainly determined 
by the DOS and not affected by the spin orbit coupling^®, 
the spin-spin correlation we study here shows rich fea¬ 
tures spatially. In general, the diagonal terms decay with 
r~'^ for n ^ 0, which is in agreement with the results in 
normal 3D metals. The spin-spin correlations are inde¬ 
pendent of the Dirac or Weyl semimetal, except for the 
spatial separation between the impurity and conduction 
electron along the z-axis. In that case, the momentum 
b/2X adds an extra phase factor to spin-spin correlation 
Jxxir) and Jyy{r), so that the spin-spin correlation is 
modified according to b. 

VI. ACKNOWLEDGEMENT 

This work is supported in part by Na¬ 
tional Basic Research Program of China 
(2014CB921201/2014CB92I203), NSEC (No.ll374256/ 
11274269), and the Fundamental Research Funds for the 
Central Universities in China. 


^ K. Novoselov, A. Geim, S. Morozov, D. Jiang, 
Y. Zhang, S. Dubonos, I. Grigorieva, and A. Firsov, 
Science 306, 666 (2004). 

^ A. H. Castro Neto, F. Guinea, N. M. R. 
Peres, K. S. Novoselov, and A. K. Geim, 
Rev. Mod. Phys. 81, 109 (2009). 

® V. N. Kotov, B. Uchoa, V. M. Pereira, F. Guinea, and 

A. H. Castro Neto, Rev. Mod. Phys. 84, 1067 (2012). 

J. E. Moore, Nature 464, 194 (2010). 

® M. Z. Hasan and C. L. Kane, 
Rev. Mod. Phys. 82, 3045 (2010). 

® X.-L. Qi and S.-C. Zhang, 

Rev. Mod. Phys. 83, 1057 (2011). 

^ Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, 
Phys. Rev. B 88, 125427 (2013). 

® Z. Liu, B. Zhou, Y. Zhang, Z. Wang, H. Weng, D. Prab- 
hakaran, S.-K. Mo, Z. Shen, Z. Fang, X. Dai, et ai, Science 
343, 864 (2014). 

® Z. Liu, J. Jiang, B. Zhou, Z. Wang, Y. Zhang, H. Weng, 
D. Prabhakaran, S.-K. Mo, H. Peng, P. Dudin, T. Kim, 
M. Hoesch, Z. Fang, X. Dai, Z. Shen, D. Feng, Z. Hussain, 
and Y. Chen, Nature Materials 13, 677 (2014). 

M. Neupane, S.-Y. Xu, R. Sankar, N. Alidoust, G. Bian, 
C. Liu, 1. Belopolski, T.-R. Chang, H.-T. Jeng, H. Lin, 
et al, Nature communications 5, 3786 (2014). 

X. Wan, A. M. Turner, A. Vishwanath, and S. Y. 
Savrasov, Phys. Rev. B 83, 205101 (2011). 

A. A. Burkov, M. D. Hook, and L. Balents, 
Phys. Rev. B 84, 235126 (2011). 


M. M. Vazifeh and M. Franz, 
Phys. Rev. Lett. Ill, 027201 (2013). 

H. Weng, C. Fang, Z. Fang, B. A. Bernevig, and X. Dai, 
Phys. Rev. X 5, 011029 (2015). 

S.-M. Huang, S.-Y. Xu, 1. Belopolski, C.-C. Lee, 

G. Chang, B. Wang, N. Alidoust, G. Bian, M. Neupane, 
C. Zhang, S. Jia, A. Bansil, H. Lin, and M. Hasan, 
Nature Communications 6, 7373 (2015). 

S.-Y. Xu, 1. Belopolski, N. Alidoust, M. Neupane, G. Bian, 
C. Zhang, R. Sankar, G. Chang, Z. Yuan, C.-C. Lee, S.-M. 
Huang, H. Zheng, J. Ma, D. Sanchez, B. Wang, A. Ban¬ 
sil, F. Chou, P. Shibayev, H. Lin, S. Jia, and M. Hasan, 
Science 349, 613 (2015). 

B. Q. Lv, H. M. Weng, B. B. Fu, X. P. Wang, 

H. Miao, J. Ma, P. Richard, X. C. Huang, L. X. Zhao, 
G. F. Chen, Z. Fang, X. Dai, T. Qian, and H. Ding, 
Phys. Rev. X 5, 031013 (2015). 

S.-Y. Xu, N. Alidoust, 1. Belopolski, Z. Yuan, G. Bian, T.- 
R. Chang, H. Zheng, V. Strocov, D. Sanchez, G. Chang, 

C. Zhang, D. Mou, Y. Wu, L. Huang, C.-C. Lee, S.- 
M. Huang, B. Wang, A. Bansil, H.-T. Jeng, T. Neu- 
pert, A. Kaminski, H. Lin, S. Jia, and M. Zahid Hasan, 
Nature Physics 11, 748 (2015). 

C. Zhang, Z. Yuan, S. Xu, Z. Lin, B. Tong, M. Z. 
Hasan, J. Wang, C. Zhang, and S. Jia, arXiv preprint 
arXiv:1502.00251 (2015). 

L. Yang, Z. Liu, Y. Sun, H. Peng, H. Yang, T. Zhang, 
B. Zhou, Y. Zhang, Y. Guo, M. Rahn, D. Prabhakaran, 
Z. Hussain, S.-K. Mo, C. Felser, B. Yan, and Y. Chen, 



Nature Physics 11, 728 (2015). 

Z. Wang, Y. Zheng, Z. Shen, Y. Zhou, X. Yang, Y. Li, 

C. Feng, and Z.-A. Xu, arXiv preprint arXiv: 1506.00924 
(2015). 

22 X. Huang, L. Zhao, Y. Long, P. Wang, D. Chen, Z. Yang, 

H. Liang, M. Xue, H. Weng, Z. Fang, X. Dai, and G. Chen, 
Phys. Rev. X 5, 031023 (2015). 

2® B. Lv, N. Xu, H. Weng, J. Ma, P. Richard, X. Huang, 

L. Zhao, G. Chen, C. Matt, F. Bisti, V. Strocov, J. Mesot, 

Z. Fang, X. Dai, T. Qian, M. Shi, and H. Ding, 
Nature Physics 11, 724 (2015). 

2"i P. W. Anderson, Physical Review 124, 41 (1961). 

2® J. Kondo, Progress of theoretical physics 32, 37 (1964). 

2® H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, 
Phys. Rev. B 21, 1003 (1980). 

21' A. Tsvelick and P. Wiegmann, Zeitschrift fiir Physik B 
Condensed Matter 54, 201 (1984). 

2® N. Andrei and C. Destri, Phys. Rev. Lett. 52, 364 (1984). 

2® F. C. Zhang and T. K. Lee, Phys. Rev. B 28, 33 (1983). 

P. Coleman, Phys. Rev. B 29, 3035 (1984). 

N. Read and D. Newns, Journal of Physics C: Solid State 
Physics 16, 3273 (1983). 

®2 Y. Kuramoto, Zeitschrift fiir Physik B Condensed Matter 53, 

O. Gunnarsson and K. Schonhammer, 
Phys. Rev. Lett. 50, 604 (1983). 


1. Affleck, Nuclear Physics B 336, 517 (1990). 

H. Ishii, Journal of Low Temperature Physics 32, 457 
(1978). 

®® V. Barzykin and 1. Affleck, Phys. Rev. B 57, 432 (1998). 
L. Borda, Phys. Rev. B 75, 041307 (2007). 

C. Gonzalez-Buxton and K. Ingersent, 
Phys. Rev. B 57, 14254 (1998). 

L. Fritz and M. Vojta, Phys. Rev. B 70, 214427 (2004). 

M. Vojta and L. Fritz, Phys. Rev. B 70, 094502 (2004). 

T. Shirakawa and S. Yunoki, 

Phys. Rev. B 90, 195109 (2014). 

“^2 A. K. Mitchell and L. Fritz, 
Phys. Rev. B 92, 121109 (2015). 

C. M. Varma and Y. Yafet, Phys. Rev. B 13, 2950 (1976). 
V. Aji, G. M. Varma, and 1. Vekhter, 

Phys. Rev. B 77, 224426 (2008). 

X.-Y. Feng, W.-Q. Chen, J.-H. Gao, Q.-H. Wang, and 
F.-G. Zhang, Phys. Rev. B 81, 235411 (2010). 

^® L. Fu and E. Berg, Phys. Rev. Lett. 105, 097001 (2010). 
K. Sengupta and G. Baskaran, 

Phys. Rev. B 77, 045417 (2008). 

H.-B. Zhuang, Q.-F. Sun, and X. Xie, 

37 pB$3^6, 58004 (2009). 

A. Principi, G. Vignale, and E. Rossi, 

Phys. Rev. B 92, 041107 (2015). 


